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Abstract. We consider a stochastic wave equation in space dimension three 
driven by a noise white in time and with an absolutely continuous correlation 
measure given by the product of a smooth function and a Riesz kernel. Let 
Pt,x(y) be the density of the law of the solution u(t,x) of such an equation 
at points (t,x) €]0,T] x M 3 . We prove that the mapping (t,x) \-> pt, x (y) 
owns the same regularity as the sample paths of the process {u(t, x), (t, x) € 
]0, T] xR 3 } established in [5]. The proof relies on Malliavin calculus and more 
explicitely, the integration by parts formula of [T7j and estimates derived 
form it. 
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1 Introduction 



In this paper we consider the stochastic wave equation in space dimension 
d = 3, 

(J^-A\u{t,x) =a(u(t,x))F(t,x) + b(u(t,x)), (1.1) 
d 

u(0, x) = v (x), x ) = 5 o(jc), 

where i G ]0, T] for a fixed T > 0, i £ R 3 , and A denotes the Laplacian on 
M 3 . The functions a and 6 are Lipschitz continuous, and the process F is the 
formal derivative of a Gaussian random field, white in time and correlated 
in space defined as follows. 
Let 

f(x)=<p(x)k /3 (x), (1.2) 

where (p is a smooth positive function and kg denotes the Riesz kernel 
kp{x) = \x\~ p , P G ]0, 2[. We shall assume that / defines a tempered measure 
and 

TTW ' ( 3) 

where = J-~ 1 (f) and J 7 denotes the Fourier transform operator. This 
condition is satisfied for instance for densities of the form (|1 .2|) with ip{x) = 
exp (— a 2 \x | 2 /2) and e]0,2[ (see [5]). 

Let P(]R 4 ) be the space of Schwartz test functions (see [15] )■ Then, on 
some probability space, there exists a Gaussian process F = (F(cp), ip G P(M 4 )) 
with mean zero and covariance functional defined by 

E[F(<p)F(i/>)) = ! ds [ dxf(x)(<p(a)*$(sj)(x), (1.4) 

where ^(s, x) = ip(s, —x). 

Riesz kernels are a class of singular correlation functions which have 
already appeared in several papers related with the stochastic heat and 
wave equations, for instance in [2], [3], [7], [8], [9]. 

We recall that the fundamental solution G(t) associated to the wave 
operator C = J^j- — A in dimension three is given by G(t) = -^(Tt, t > 0, 
where cr^ denotes the uniform surface measure on the sphere of radius t G 
[0, T], hence with total mass 4irt 2 . 
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The properties of G(t) together with the particular form of the covariance 
of the noise play a crucial role in giving a rigourous formulation to the initial 
value problem (jl.ip . 

Here, we shall follow the same formulation as in [5] which for the purpose 
of existence and uniqueness of solution of (jl.ip introduces a localization of 
the SPDE by means of a set related with the past light cone, as follows. Let 
D be a bounded domain in IR 3 . Set 

K^(t) = {yeR i :d(y,B)<a(T-t)}, t € [0,7], (1.5) 

where a > 1 and d denotes the Euclidean distance. Then, a solution to the 
SPDE (jl.ip in D is an adapted, mean-square continuous stochastic process 

u(t)l K o^, t € [0,T]J with values in L 2 (M 3 ), satisfying 

<t, 0lx5>(t)(0 = ljcg(*)(0 (^ G W * + G(t) * £ ) (•) 

+ 1 A-D (t )(-) / / G(t-a l --y)o-(«(s,y))l x n W (y)Af(ds 1 dy) 

+ 1k?w(')J dsG ^- s ">* ■))!*?«(•)) . (1-6) 

a.s., for any i € [0, T], where we consider the stochastic integral defined in 
[I] and M denotes the martingale measure derived from F (see [2]). 

The following result is a quotation of Theorem 4.11 in [5] and will be 
invoked repeatedly in this paper. 

Theorem 1.1 Assume that: 

(a) the covariance density is of the form with the covariance factory 
bounded and positive, cp € C 1 (M. d ) and Vip € C^(M d ), for some 6 €]0, 1]; 



(b) the initial values vq, vq o^e such that vo € C 2 (M 3 ) ; and Avq and vq are 
Holder continuous with orders 71,72 S]0, 1], respectively; 

(c) the coefficients a and b are Lipschitz. 

7hen for any q E [2, 00 [ and a € ]0, 71 A 72 A A there is C > such 
that for (t,x),(i,y) e[0,T] xD, 

£ ([«(*, z) - u(t, y)| 9 ) <C(\t- i\ + \x- y\) aq . (1.7) 
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In particular, a.s., the stochastic process (u(t,x), (t,x) € [0, T] x D) solu- 
tion of l\1.6\) has a-Holder continuous sample paths, jointly in (t,x), and 



sup E(\u{t,x)\ q ) < oo (li 
(<,i)e[o,T]xD 



for any q G [1, oo[. 



In this paper we are interested in studying the properties of the density 
of the solution of (|1 .6j) as a function of (t, x) g]0, T] x D, where D is a 
bounded subset of M 3 . We shall denote this density by Pt, x (y)- We shall 
prove that (t, x) — > Pt,x(y) is jointly Holder continuous, uniformly in y on 
compact sets. 

This question is trivial in the very particular case where the initial con- 
ditions Vq, vq and the coefficient b vanish, and the coefficient a is a constant 
function. In fact, with these assumptions and a = 1 the solution to the equa- 
tion (|l,6p is a Gaussian process, centered, stationary in the space variable, 
and with 

rt r sin 2 ( S |e|) 



of := E\u(t,x)\ 2 = [ ds [ (j,(d£) 
Jo Jr 3 



From the expression pt,x(y) = 1 1 exp ( — ipr) it is not difficult to prove 
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that 

\pt,x(y) -Pt,x(y)\ ^ c{t ,ti,B) \t-t\, 

for < t < t < i < h, D C M 3 . 

However, in the general situation that we are considering in this article, 
the problem becomes much more involved. 

Suppose that vq, do are null functions, assume also that the covariance 
of the process F is given by (|1.4j) with dxf(x) replaced by T(dx), where 
T is a non-negative, tempered, non-negative definite measure. Set /x = 
^ r_1 (r). We introduce an assumption on fi, denoted by (H^), saying that 
J R3 ^^gp^ < oo, for some value of r/ €]0, 1]. Assume that the coefficients 
a and b are of class C 1 with bounded derivatives and that (H^) holds for 
some 7] E (0, 1). Then, the existence of the density pt, x at any fixed point 
(t, x) e]0, T] xD has been established in [13]. Moreover, assuming that a and 
b are C°° functions with bounded derivatives of order greater or equal to one, 
and that (H^) holds for some r\ G (0, |), it is proved in [TJ] that y ^ Pt,x{y) 
is a C°° function. We refer the reader to [16J for results on applications of 
Malliavin calculus to the analysis of probability laws of SPDEs. 

In [12], it is shown that the extension of Walsh's integral introduced in 
[3] does not require for the integrands any stationary property in the spatial 
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variable. As a consequence of this fact, the results of [3], [13], [2] and [16] 
concerning the stochastic wave equation can be formulated with non null 
deterministic initial conditions. In addition, the solution of the equation 
in this setting coincides with the solution to (jl.6p . Furthermore, in the 
particular case of absolutely continuous covariance measures T(dx) = f(x)dx 
satisfying (|1.3|) the existence and smoothness of the density pt^ x are proved 
in [12] under the weaker assumption (Hi). 

Hence, on the basis of the above mentioned references and remarks, we 
can write the next statement, which together with Theorem 11.11 are the 
starting point of our work. 

Theorem 1.2 Assume assumptions (a) and (b) of Theorem li.il Suppose 
also that a and b are C°° functions with bounded derivatives of order greater 
or equal to one, and inf{|c(z)|, z € M} > gq > 0. Then, for any fixed 
(t,x) g]0, T] x D, the law of the real valued random variable u(t,x) solution 
to M.6\) has a density pt >x € C°° . 

The main purpose of this paper is to prove that with the assumptions of 
this theorem, for any y € M, the mapping 

(t,x) E]0,T] xR 3 ^p t:X (y) 

is jointly a- Holder continuous with a e]0, inf (71,72, (see Theorem 

O in Section ED. 

For stochastic differential equations and some finite-dimensional stochas- 
tic evolution systems with an underlying semigroup structure one can find 
results of this type for instance in ([6]). For SPDEs the problem has not 
been yet very much explored. To the best of our knowledge, this issue has 
only been studied for the stochastic heat equation in spatial dimension d = 1 
in [10] and for the wave equation with d = 2 in [9] (see [JJ and [7] for the 
existence and regularity of the density for these two types of SPDEs). It is 
worthy noticing that in these two references, the Holder degree regularity of 
Pt,x(y) in (t,x) is better than for the sample paths of the solution process 
u(t, x), while in the equation under consideration we obtain the same order. 
As it will become clear from the proof, the reason is the rather degenerate 
character of the fundamental solution of the wave equation in dimension 
three. 

The method of our proof is based on the integration by parts formula 
of Malliavin calculus, as in the above mentioned references. We next give 
the main ideas and steps of the proof. Fix y € R and let (g n ,y,n > 1) 
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be a sequence of smooth functions converging pointwise to the Dirac delta 
function 5s y \. Fix t,i e]0,T], x,x € D and assume that we can prove 

sup \E\g n} y{u(t,x))-g n , y {u(t,x))}\<c\\t-t\^ + \x-xf 2 ], (1.9) 

for some Pi,^2 > 0, where K C R. Then, since Pt,x(y) = E [S{ y } (u(t,x))~\ 
(see Theorem 1.12 in [T7| for a rigorous meaning of this identity), by passing 
to the limit as n — > oo, we will have joint Holder continuity of the mapping 
(t,;c) g]0, T] x D i— > pt,x(y) € R with degree /3\ in i and (3% in x, uniformly 
in y G X. 

An estimate like (jl.9p is obtained by the following procedure. For sim- 
plicity we write g instead of g n ,y We first consider a Taylor expansion of 
g(u(t,x)) around u(t,x) up to a certain order tq chosen in such a way to 
obtain optimal values of j3\ and fy. Then for any r < ro, we estimate terms 
of the type 



E 



(u(t,x))(u(t,x)-u(t,x)y 



and the term corresponding to the rest in the Taylor expansion, whose struc- 
ture is similar. For this we use the version of the integration by parts formula 
for one-dimensional random variables given in Lemma 2, page 54 of |17j (see 
also Equations (2.29)-(2.31) of [TT]) which we now quote as a lemma. 

Lemma 1.3 On an abstract Wiener space (f2,H,P), we consider two real- 
valued random variables £ and Z such that £ G B°° ; H-OCIIh 1 ^ ^p>2L p (Q), 
Z € B°°. Let g be a function in C , for some r > 1. Denote by g the 
antiderivative of g. Then, the following formula holds: 

E (gV(Z)z) =E(g(0H r+ i(Z,0) , (1.10) 

where H r , r > 1, is defined recursively by 

ZDi \ 

'r 

ZDi 



n 

H r+1 (Z,0 = 5 (^H r (Z,0j^rJ ,r > 1. 

In this Lemma, 5 stands for the adjoint operator of the Malliavin derivative, 
also termed divergence operator or Skorohod integral and we have used the 
notations of [11] and [16], as we shall do throughout the paper when referring 
to notions and results of Malliavin calculus. 
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The abstract Wiener space that we shall consider here is the one asso- 
ciated with the Gaussian process F restricted to the time interval [0, T], 
as is described in [16] , Section 6.1. For the sake of completeness and its 
further use, we recall that H := H T , U T = L 2 ([0, T];H) and that H is the 
completion of the inner product space consisting of test functions endowed 
with the inner product 



n 



dxf(x) ( cp * if) ) (x). 



Assume that the function g in Lemma 11.31 is bounded. From (jl.lOp it 
clearly follows that 



E 



9 {T) m] < \\~g\U\H r+1 (Z,0hi 



(«)■ 



1.11 



Furthermore, as a consequence of the continuity property of the Skorohod 
integral and the assumptions on £, for any r > 1, k > 1 and p € (1, oo), 

\\H r (Z, C)\\k,p < C||-^||fc+r,4 r p 

(see Corollary 4.1 in [10J). Consequently, under the previous assumptions 
from (jl.lip we obtain 



E 



^ C||5l|oo||-^||r+l. 



4r-+l ■ 



(1.12) 



Let us recall that for a natural number k and a real number p £ [1, co[, 

k 



\k,p 



\z\ 



r=l 



We shall apply (|1 . 12|) mainly to £ := u(t, x) and Z := (u(t, x) — uii, x)) r , for 
natural values of r. Under the hypotheses of Theorem 12.11 the assumptions 
of Lemma 11.31 are satisfied (see [13] and Chapters 7 and 8 of [16] ) . Thus we 
face the problem of giving upper bounds for || (u(t,x) — u(i,x)) r Kr+i^+i. 

Malliavin derivatives of the solution of (jl.ip satisfy evolution equations 
(see |16| . Theorem 7.1 and |14j). Indeed, for x € D, and a natural number 
k > 1, (D k ^u{t,x),(t,x) G [0,T] x D) is a ftf^-valued process satisfying 
D T ^u(t, x) = for r > t, and for r < t it is the solution of the evolution 
equation 



D k r Mt,x) = Z k (t,x) 
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+ / / G(t-s,x-y) 
Jo Jr3 

x (r k (a,u(s,y)) +a > (u(s,y))D^u(s, y)) M(ds,dy) 

+ f ds [ G(s,dy)(T k (b,u(s,y)) + b'{u(t-s,x-y))D^u{t-s,x-y) 
Jo Jr 3 v 

(1.13) 

In this equation, (Z k (t, x), (t, x) E [0,T] x D) is a Wf! fc -valued stochastic 
process and for a given function g € C k and a random variable X G B fc ' 2 , 
r fc ( 5 , X) = D k (g{X)) - g'(X)D k X. 
The solution of (|1 . 13j) satisfies 

sup \\u{s,y)\\ k < +00 (1.14) 
(t,x)e[0,T]xD 

for any p € [1, oof (see Theorem 7.1 in [16 J. 

In the next section, we shall make use of the explicit form of (|1.13p for 
k = 1. In this case 

Z^(t, x) =G(t-r,x- *)a(u(r, *)) (1.15) 

and T^X) = 0. 

With some effort, using the tools on stochastic integration of Hilbert- 
valued processes developed in [16] it can be proved that the conclusions 
of Theorem II .11 also apply to the 7^® fc -valued stochastic process solution to 
(fLT3j) . More precisely, for any k > 1, q G [2, oof and a G]0, 71 A7 2 A^A 1 ±^[, 
there is C > such that for (t, x), (i, y) G [0, T] x D, 

\\u(t,x)-u{i,y)\\ k>q < C(\t-t\ + \x-y\) a . (1.16) 

Hence, with the Holder continuity property on u(t, x) and its Malliavin 
derivatives we may be able to prove (|1.9|) for specific values of (3\ , ■ 

We shall fix what is the top order r$ in the Taylor expansion of g{u(i.x)). 
Clearly, the lower exponents (3i should come from the first order term. How- 
ever, in the examples studied so far, terms of first and second order give the 
same exponent. For the equation (II. 6p the situation is different. Already 
at the first order level of the expansion, we shall see that the contribution 
of the pathwise integral involving the coefficient b is of the same order than 
the Holder continuity exponent given in Theorem 11.11 Clearly, the second 
order term would provide twice the Holder continuity degree. Therefore, a 
Taylor expansion of first order gives the best possible result. However, to 
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conclude whether the regularity of the density pt,x in (t, x) is the same as 
that of the sample paths of u(t,x), we have to check that the contribution 
to the first order term in the Taylor expansion of the stochastic integral is 
not worse than that of the pathwise integral. This explains the strategy of 
the proof of the main result in the next section. 



2 Main Result 

Throughout this section D denotes a fixed bounded domain of R 3 and C will 
be any positive finite constant. We assume that (|1.3|) holds. Our purpose is 
to prove the following theorem. 

Theorem 2.1 Assume that: 

(a) the covariance density is of the form il.2\) and the covariance density 
factor ip is bounded and positive, ip G C 1 (]R c( ) and V</? G C^(M. d ) for 
some 5 G]0, 1]; 

(b) the initial values vq, vq are such that vq € C 2 (R 3 ), and A-uo and vq 
are Holder continuous with orders 71,72 G]0, 1], respectively; 

(c) the coefficients a and b are C°° functions with bounded derivatives 
of order greater or equal to one, and there exist o~q > such that 
mf{\a{z)\,z £l}> a . 

Then the mapping 

(t,x) G]0,T] x D 1— > pt,x(y) 

is a-Holder continuous jointly in (t,x) with a g]0, inf(7i, 72, ^S^, uni- 
formly in y G R. 

Proof: Fix y G R and let (g n ^ y ,n > 1) be a sequence of regular functions 
converging pointwise to 5^ y y as n — > 00; for example, a sequence of Gaussian 
kernels with mean y and variances converging to zero. We may assume 
that the corresponding antiderivatives g niV are uniformly bounded by 1. To 
simplify the notation, we shall write g instead of g n ,y 

1. Time increments 

For (t, x) G [0, T] x D we consider the Taylor expansion 

E [g (u(t + h, x)) - g (u(t, x))] = E [g' (u(t, x)) (u(t + h,x)- u(t, x))] 



+ E 



g" (u(t,x,h)) (u(t + h,x) — u(t,x)) , (2.1 



S 



where h > and u(t, x, h) denotes a random variable lying on the segment 
determined by u(t + h, x) and u(t, x). 
First order term 
Set 

Ti(t, x,h) = \E [g' (u(t, x)) (u(t + h,x)- u(t, x))] \ . 
We aim to prove that 



sup T x (t,x,h) < Ch a , 

(t,x)e[0,T]xD 



(2.2) 



2-/3 A 1+6 \ 



with a G]0, 71 A 72 A 2 , , 2 

By using equation (jl.6p . we write T±(t, x, h) < 5^t=l ^~i,«(*> x ' w hh 



</ («(t, x)) ( — G(t + h)*v + G(t + h)*v ) (x) 



G(t) * v + G(t) * v (x) 



t+h 



T\p{t, x, h) = E 1 g' (u(t,x)) / ds / G(i + h — s, dz)b(u(s, x — z)) 



T 1>3 (t,x,h) 



E 



ds I G(t — s, dz)b(u(s, x — z)) 

g (u(t,x)) 

t+h 



G(t + h — s, x — z)a(u(s, x — z))M(ds, dz) 
G(t — s, x — z)a(u(s, x — z))M(ds, dz) 



In fact, by our choice of (t,x) all the indicator functions in (|1.6p take the 
value 1. 

We shall apply repeatedly the inequality (|1.12p with r = 1, £ := u(t,x) 
and different choices of Z. To begin with, we take 



Z := \-G(t + h)*v + G(t + h) * v j (x)-l—G(t) * v + G(t) * v J (x). 

Since Z is deterministic, ||Z||fc,p = \Z\, for any k and p. Then, applying 
(|1.12j) and Lemma 4.9 in [5] yields 

sup Ti,i(t,x,h) < Ch? 1 ^ 2 . (2.3) 
(i,i)e]o,r]xD 
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We next study the term T\^{t, x, h). Let 

rt+h 



Ti t 2,i(t,x,h) = E g' (u(t,x)) / ds G(t + h — s , dz)b(u(s , x — z)) 
L Jt Jr 3 - 1 

We apply (|1 . 12j) to Z := J* + ds L 3 G(t+h — s,dz)b(u(s,x — z)) and consider 

the measure on [t, t + h] x M 3 given by dsG(t + h — s, dz) with total mass 4p 
and an arbitrary p£ [1, oof. By applying Minkowski's inequality, we obtain 



t+h 



< 



ds I G(t + h — s,dz)b(u(s,x — z)) 
t+h 



2-P 



ds / G(t + h — s, dz) \\b(u(s, x — z))\\ 2 „ 
Jk. 3 

By the chain rule of Malliavin calculus, 

\\ b (u(s,y))\\ 2)P < C(l + \\u(s,y)\\ 2 ^ p + \\u(s,y)\\l >2p 
Consequently, 

sup \\b(u(s,y))\\ 2 „ < oo. 

(s,y)e[0,T]xD 2T 



(2.4) 



where for a bounded set D C M 3 and a > 0, we denote by D a = {z G 
R 3 ;d(z,D) < a} and we have applied (|1.14p . 
Thus, we have proved 



sup 

(t,i)£]0,T]xD 



t+h 



ds / G(t + h — s,dz)b(u(s,x — z)) 



<Ch 2 



2-P 



for any p E [1, oof, and consequently, 



sup r 1)2 ,i(i,x,^) < C/i 2 . 

((,i)e]0,T]xD 



(2.5) 



Set 

Tl,2,2(t,X,k) 



E 



g' (u(t,x)) ds b(u(s,x—z))[G(t+h—s,dz)—G(t—s,dz)] 



that according to (| 1 . 12[) we can bound as follows, 

Ti, 2 ,2(t,x,h) <C ( ds I b(u(s,x-z))[G(t + h-s,dz)-G{t-s,dz)] 
Jo Jk 3 



2,4 2 
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We have 



Bi(0) 



Bi(0) 



ds / b(u(s, x — z))G(t + /i — s, dz) = J ds(t + h — s) 
G(l, dz)b(u(s, x — (t + h — s)z)), 

ils j b(u(s, x — z))G{t — s, dz) = I ds(t — s) 

G(l,dz)b(u(s,x - (t- s)z)), (2.6) 
as can be easily checked by applying the change of variables z t— > t+ ^_ s and 
P G [l,°o[ 

ds [ b(n(s.x - z))[G(t + h- s,dz) - G(t- s,d;,)\ 
G(l, dz)b (u(s, x- (t + h- s)z) 



z i — > r^j, respectively. Then, by the triangular inequality we obtain for any 



2.p 



< Ch 



+ C 



ds 

^Bi(O) 



2,p 



ds 

^Bi(O) 



G(l, dz) [6 (u(s, cc — (t + /i — s)z) — 6 (u(s, cc — (t — s)z)] 



2,p 



(2.7) 



For the study of the first term in the right-hand side of (|2.7p . we apply 
Minkowski's inequality and then (|2.4p . This yields 



Bi(0) 



G(l, dz)6 (it(s, a; - (t + h - s)z) 



< C. 



2.p 



The Lipschitz property of b and (TOD, ([116]) . ([T7H]) yield 

||6(it(s,ac- {t + h- s)z) - b(u(s,x - (t - s)z)\\ 2 < Ch a , 
with q g]0, 7i A72A^^A j ^^[. Consequently, after having applied Minkowski's 



inequality we see that the second term of the right-hand side of (|2.7p is 
bounded by Ch a , uniformly in (t,x) € [0, T] x D. 
Thus, we have proved 

sup T 1)2j2 (t,x,h) < Ch a , 

(t,i)6[0,T]xD 
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and along with (|2.5p we obtain 



sup T 1>2 {t,x,h) <Ch° 

(t,x)e[0,T]xD 



(2. 



with a G]0,7i A 72 A ^ A ±±^[. 

Let us remark that in [10] and [9] the contribution of the term analogous 
to Ti^(t, x, h) is a power of h of higher order than the Holder degree of 
the solution. In fact, for the heat equation and the wave equation in spatial 
dimension two, by integration of the increments G(t + h — s, dz) — G(t — s, dz) 
we get powers of h. For the wave equation in dimension three, such an 
approach is not possible. Instead, "increments" of G(t — s,dz) are transfered 
to increments of b(u(s,x — z)) (this is the role played by the change of 
variables that we have performed to obtain (|2.6h ) and after this, we can 
conclude by applying the Lipschitz property of b and the Holder continuity 
of the sample paths. 

The inequality (|2.8p tell us that we are not going to improve the Holder 
degree of the mapping t e]0, T] h- > pt iX (y) hi more than the given a. But 
it might happen that the contribution of Tis(t, x, h) makes the overall esti- 
mation worse. We next carry out a careful analysis of this term and prove 
that its contribution in terms of powers of h is the same as T\^it, x, h). 

We write Ti^(t, x, h) = 11,3,1 (i, x, h) + T\^^{t, x, h) with 



11,3,1 (*, a, h) 

rt+h 



E 



g'(u(t,x)) 



a(u(s, y))G(t + h — s, x — y)M(ds, dy) 



T li32 (t,x,h) 



E 



g'(u{t,x)) 



II a(u(s,y))[G(t + h - s,x - y) - G(t - s,x - y)] M(ds,dy) 
/o Js. 3 

The term Txsi(t,x,h) vanishes, since the random variable g'(u(t,x)) is 
adapted to the natural filtration generated by the martingale measure M. 

In contrast with Ti,2, for the analysis of % (t, x, h) we do not start 
by applying (11.121) . which actually would lead to worse results (see Remark 
2.2|) ; instead, we apply Proposition 3.9 of [16]. Since the mathematical 
expectation of a Skorohod integral is zero, we obtain 



Ti,3,2(t,x,h) = E^g"(u(t,x))D.^u(t,x),a(u(.,*)) 

x [ G (t + h _ x _ *) _ G {t _ x _ *)] i ]M (.)^ 
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7~Lt 



E \g"{u(t, x))(D. j# u(t, x),a(u(., *)) 
x [G(t + h-.,x-*)-G(t-.,x-*)} l m (•)) 

For any (t, x) G [0, T] xR 3 , we define 

B h > ,(i,x) = a(u{.,*)) [G(t + h-.,x-*)-G(t -.,£-*)] , 
With this notation, and by applying (|1.12p to T\ 3 2 (t, x, h) we see that 



T li3i2 (t,x,h) < C 



D.^u(t,x),B h .,*(t,x] 



1~Lt 



3,4 3 



(2.9) 



We shall consider norms || ■ ||fc iP for arbitrary k > 1 and p£ [l?oo[, instead 
of || • H343. By virtue of (|1.13p and (|1.15p . we write 



D.^u(t,x),B h .,*(t,x) 



k.p 



< C (\\B^(t,x;i,x)\\ kjP + \\B%(t,x;i,x)\\ kiP + \\B^{t,x;i,x)\\ kj ,) , 



where 



B^(t,x;i,x) = {G(t - .,x- *)<r(u(., *)), B h ., x) 



T~Ct 



B$(t,x;t,x) = ( f f G(t-s,x-y)a'(u(s,y))D^u(s,y)M(ds,dy), 
x Jo Jrz 

B h .,*(t,x)) , 

B% (t, x; t, x) = I I ds I G(s,dy)b' (u(t — s,x — y))D #u(t — s,x — y), 
\ Jo Jr 3 



B h .,*(t,x] 



Consider the change of variables (s,y) — > yt — s,-^); Fubini's theorem 
along with Minkowski's inequality yield 



B$(t,x;t,x] 



k.p 



ds(t-s) / G(l,dy)b'(u(s,x - (t-s)y)) 



x (D ^u(s,x — (t — s)y), B h .,*(t,x) 



k.p 
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< C\W 



ds / G(l,dy) 



D.*u(s, x - (t - s)y),B h ., *(t, x) 



Hj 



k-p 



< CI ds sup 

to yeK°(s) 



D,u{s,y),B h .,*(i,x t 



Hj 



k,p 



The last inequality is obtained as follows. By definition of the sets K^(t), 
it is obvious that for any (t,x) € [0,T] x D, x belongs to K®{t). Since the 
support of the measure G(s, dy) is the boundary of the ball centered at zero 
and with radius s, the y-variable in the above integrals belongs to K®{s). 
Hence, by GronwalFs lemma 



sup 



D±u(t,x),B h Jt,x) 



k,p 



<C sup [\\B?(t,x;t,x)\\ k , p + \\B%(t,x;t,x)\\ kiP 



Our aim is to prove that 



sup 

(t,i)£[0,T]xD 



D j# u(t,x),B h Jt,x) 



< Ch° 



(2.10) 



(2.11) 



k.p 



with a G]0,7i A 72 A A ^p[ 3 for k = 3. This is done recursively on 
k = 0, 1, 2, 3, where by convention || • ||o, p = || • 

To illustrate the method and simplify the presentation, we shall consider 
in (|2.10p the norm || • instead of || • ||3 P . That is, we shall deal only with 
derivatives up to the first order. Thus let as first prove (12. lip for k = 0, that 
means for the L p (Q)-norm. For this, we start by studying the L p (f2)-norm 
of B%(t,x;t,x). 

To shorten the notation, set T>(s,y;t,x) = (D^u(s,y), B\{f,x))% T . 
From Burkholder's inequality it follows that 

\\B%(t,x;t,xW p 



= E 
< CE 



JM 3 
t 



G(t - s, x - y)a' (u(s, y)) V(s, y; t, x)M(ds, dy) 



ds 



dy / dyG(t - s, x - y)a (u{s, y)) f{y - y) 



xG(t-s,x- y)a (u(s, y)) V(s, y; t, x)V(s, y; t, x) 



< C / ds sup 

/o yeK°(s) 



E 



D.,Ms,y),B h Jt,x) 



(2.12) 
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Therefore, (|2.10p with k = and Gronwall's lemma yields 



sup 



D^u{t,x),B h ^{t,; 



Hn 



< C sup 



B^t,x;t,x) 



with a constant C independent of t and x. Hence, we can fix t = t and x = x 
in the preceding inequality and obtain 



sup sup 

te[o,T] xeiff(i) 



D^u(t,x),B h M,x) 



< C sup sup 

p te[o,T] xexfW 



where B^(t,x) stands for B±(t,x;t,x). 

By the very definition of the inner product in TLt we have 



(2.13) 



B\(t,x) 



E 



dr dU d V G(t-r,x-ti)a(u(r,0)f(ti-ri) 



X [G(t + h — r, x — rj) — G(t — r, x — 77)] a(u(r, 77)) 



We consider each one of the terms in the difference of the right-hand side of 
this inequality and apply respectively the change of variables 



(£,77) -> ( x - f, (x - 77) 



t — r 
t + h — r 



, (£,77) -> (x - £,x - 77). 



With this, the increments in time of the measure G are transfered to incre- 
ments of a and /. More precisely, we obtain 



B%(t,x) < C{T 1 . 3:2tl {t,x,h)+T 1:3 . 2t2 (t,x,h)) 
v 



with 



2l j3A i (*, x,h) = E ! dr I ! G(t-r,d£)G(t-r, 
JO Jr 3 Jr 3 



dr])a(u(r, x — £)) 



t + h-r (t + h-r 
x — / — : V ~ £ 



t — r 



t - r 



a(u(r, x — 77)) — a(u(r, x — 77)) 



t - r 



T x ^ 2 {t,x,h)=E f dr f I G(t-r,dt)G(t 
Jo Jr 3 Jr 3 



r, dr])a(u(r, x — £)) 



x a(u(r, x — rj)) 



t + h — r „ ( t + h — r 



t — r 



t — r 



v-a-f(v-o 
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For the analysis of ^1^2,1 (t, x, h) we consider the measure with support on 
[0,T] x B t _ r (0) x £ t _' r '(6) defined by 



t + h-r , (t + h-r 

f — : »7-£ 



v(dr;d£,dr]) := dr G(t - r,d£)G(t - r,drj) 

t — r \ t — r 
Following the steps of the proof of Lemma 6.3 in [5] we obtain 

sup / / / u{dr;d^,drj) < 00. 

0<t<t+h<TJ0 iB t _ r (0) iBi- r (0) 

Then, we can write 

ct 



(2.14) 



T 1>3j2>1 (t,x,h) = E [ [ [ 

JO JB t -r(0) Jh 



u(dr; drj, d£)a(u(r, x — £)) 



(0) JB t - r {0) 

a(u(r, x 77)) — a(u(r, x — rj)) 



and apply Holder's inequality with respect to the measure v{dr; dr], d£). This 
yields 



Tl, 3 ,2,l(t,X,h) <C I I I 
JO JB t - r (0) Jh 

a(u(r,x - £)) 



u{dr; dr], d£) 



x E 



t-r(0) JB t -r(0) 

a(u(r, x — rj)) — a(u(r, x — ri)) 

t — r 



We now apply Schwarz' inequality to the factor containing the expectation. 
Since the coefficient a is a Lipschitz function, by (jl.7p and (jl.8p we obtain 



Ti^2,i{t,x,h) <C sup E 
H=t-r V 



a{u(r, x — rj)) — a(u(r, x — rj)) 



t - r 



2p\ 2 



< C sup 

\r)\=t-r 



h 



t - r 



< Ch ap , 



(2.15) 



with a s]0, 71 A 72 A A and a constant C not depending on (t, x) G 
[0, T] x D for an arbitrary D. 

To study ^1^2,2 (t,x,h) we consider the measure on [0,t] x B t - r (0) x 
J B t _ r (0) given by 



*i h (dr; d£, dr]) = dr G(t - r, d£)G(t - r, dr]) 
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t - r 



t — r 



We also consider two additional measures with the same support as ]i h (dr; d£, drj) 
obtained by applying the triangular inequality to the expression 



t + h — r „ / 1 + h — r 



t — r 



-f 



t — r 



They are defined by 



rf(dr; d£, dr.) = drG(t - r, £)G(t - r, dr,)- f( V - £), 

t — r 



\i\ (dr; d£, dr]) = drG(t - r, £)G(t - r, dr]) 



t + h- r 
t — r 



(2.16) 



With these new ingredients, Ti 3,2,2(t, x,h) < C (Ti 3 2 ,2,1 + T\ z ,2,2.2) , where 



^1,3,2,2,1 — E 



1,3,2,2,2 



E 



t 

JB t -r(0) JB t -r(0) 
t 



Hi(dr; d£, dr])a(u(r, x — £))a(u(r, x — rj)) 



JBt-r(0) JB t -r(0) 



y!\{dr\ d£, drj)cr(u(r, x — £))a(u(r, x — rj)) 



We next check that f Q J Bt (0) Je t (0) ^1 (^ r ' ^v) < Indeed, ow- 
ing to (jl.2l) and by the change of variable r 1— > t — r, we have 



l4(dr;d£,drj) 



JB t -r(p) JB t -r(0) 
1 dr 

<Ch ■ 
= Ch 



r JB t _ r (0) JB t - r (0) 
1 dr 



G(r,d^)G(r,dr,)k^- V ) 







|^G(0(Q| 2 

d4 j^Z^ < C/l, 



uniformly in t € [0, T], where in the last inequality we have applied Lemma 
2.3 of [5J with 6=1. 

Consequently, Holder's inequality, the linear growth of the coefficient a 
and the property (|1.8|) yield 



T X , 3,2,2,1 < ChP, 



(2.17) 
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uniformly in (t, x) £ [0, T] x D. 

The next step consists of proving that (1% {dr\ d£, drj) defines a finite mea- 
sure as well, and in giving an estimate of its total mass in terms of powers 
of h. For this, we consider the inequality 



t + h-r 



t - r 
x kp{r] - 

t + h — r 



v-a-f(v-o 



< 



t + h-r 
t - r 



rj-n -<p(v~0 



+ 



t — r 



t + h-r 
t — r 



(2.18) 



which is a consequence of (|1.2p and the triangular inequality. The properties 
of (p together with Lemma 2.3 in [5] yield 



dr 



o 



X if 

< Ch 

< Ch 



Bt-r(0) JBt-r(P) 

t + h-r 



G(t-r,d^)G(t-r,drj) 



t + h-r 



t - r 



t — r 
1 dr 
t - r 



v - £ - <p(v - 



M 7 ? - 



B t _ r (0) JB t - r (0) 



G(t-r,dt)G(t-r,dri)k(,(Ti-Z) 



(2.19) 



uniformly in t € [0, T] . 

Let us now consider the contribution to ^(dr; <i£, dr]) of the second term 
of the right-hand side of (12TT8D . Let a e]0, 1[, G]0, 2{ with a + (5 e]0, 2[. 
By applying Lemma 2.6 (a) of [5] with b := a, a := 3 — (a + /?), c := /i, 



u := rj — x := t^-, we obtain 



G(t-r,d£)G(t-r,drj) 



B t -r(0) -/Bt-r(O) 

t + h — r 



t + h-r 
t - r 



t + h-r 



-0 



t - r 



v - C J - M 7 ? - 6 



< IM 



o 



* , t + h-r 

dr 

i — r 



B t _ r (0) ^Bt-r(O) 



G(t-r,d^)G(t-r,dri)Dk ( r? - 







t — r 

x / dwk &+/3 (r] - £ - /iu>) 



S t _ r (0) ^St-r(O) 



G(t - r, d£)G(t - r, dr?) 



t — r 



(2.20) 
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where we have set Dg(x, y) := g(x + y) — g{x) for a function g : M 3 — > M. 

Our next purpose is to prove that the last integral in the above expression 
is bounded, uniformly in t, t + h £ [0, T]. For this, as in Lemma 6.4 of [5] we 
split the integral on the w- variable in the last expression into the sum of two 
integrals: on a finite ball containing the origin and on the complementary 
of this set. In this way we obtain as an upper bound of 



* , t + h-r 

dr 

t — r 



B t - r (0) JB t - r (p) 



G(t - r, dQG(t - r, dr,) 



x / dwka+piv ~ £ - hw) 



Dk 3 - & w, 



t - r 



the sum of the three terms: 



(i) 



* , t+h-r 

dr 

t — r 



Bt- T (0) JB t -r{0) 



G(t - r, dQG(t - r, dr,) 



dwka+piv ~ £ - hw)k 3 - & (^w + 



V 



r(2) 







B 2 (0) 

* , t+h-r 

dr 

t — r 



t - r / 

G{t - r, d£)G(t - r, dr,) 



B t _ r (0) JB t -r(0) 



r(3) 



dwka+piv - £ - hw)k 3 _ & (w) 

B 2 (0) 

1 , t + h-r 
dr 



o 



t — r 



St-r(0) JBt-r{0) 



G(t - r, d£)G(t - r, dr,) 



B 2 (0) c 



dwka+pir, — £ — hw) 



Dkz-a. w, 



t — r 



Consider the change of variable w w + -fl- and then r, t— » *~t"^ ' r, that 



we apply to By Fubini's theorem we obtain 



B 3 (0) 



t + h-r 



t — r 



1 , t + h-r 

dr 

t — r 

r, — £ — hw 



o 



B t _ r (0) JB t _ r (0) 



G(t-r,d£)G(t-r,dfj) 



[ dwk 3 - a (w)[ drj [ G{t-r,d£)G{t-r + h,dr,) 

JBz(O) JO JR 3 JR 3 



dw 

'B 3 (0) 

X ka+piv - £ - hw). 
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The properties of the Fourier transform and the expression of this oper- 
ator applied to Riesz kernels yield, after regularization of G, 

/ / G(t-r,dZ)G(t-r + h,drj)ka+p(v-t-hw) 

JR 3 J R 3 

= [ TG[t - r)(t)FG(t - r + /0(|)*3_ (a+/3) (. - W)(£)- 

JR 3 

Hence by applying Schwarz's inequality, the last integral is bounded by 

\FG(t)(Z)\ 2 k 3 _ (a+l3) dZ, 



sup 

te[o,T] Jr 3 

which is known to be finite whenever a+(3 £ [0, 2] (see for instance Equation 
(2.5) in 0). 

Since k%-a{w) is integrable in a neighbourhood of the origin for any 
a > 0, we finally obtain is bounded uniformly in t, h € [0, T]. 

(2) 

Similar but simpler arguments show that the same property hold for J| . 



For any A € [0, 1] set tjj(X) = k^a yw + A^J. It is easy to check that 
l^'MI < Cki-a [w + X-j^-p )• Moreover, for \w\ > 2, and \rj\ = t — r, 



tu + A- 



t - r 



> 



t — r 



>H-i>V, 



by the triangular inequality 
Thus, 



if } <C I dr 



t + h-r 
t — r 



B t _ r (0) ^B t _ r (0) 



G{t-r,d£)G(t-r,drj) 



\ ka+piv hw)dw / /c 4 _ 5 ( u; + A—^— ) 

J(B 2 (0))° Jo V 



f/A 



< C 



(B2(0)) e 



dwki-a(w) 



(2.21) 



For d €]0, 1[ the integral J( B2 ( ))c ^ ^4-a(^) is finite. Moreover, if d + /3 G 
]0, 2[ the last integral in (|2.21[) is also finite, owing to Lemma 2.3 in [5] 

(3) 

applied to the value 6 = 1. This lead us to conclude that I{ is bounded 
uniformly in t, h £ [0, T]. 
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Summarizing, as a consequence of (|2.19p . (|2.20|) and the preceding dis- 
cussion, we have proved that 



sup 

te[o,t] Jo JB t -r(a) JB t -r(o) 



(i%(dr;d£,dr)) < Ch & , 



(2.22) 



with a G]0, ^£[. 

We can now apply Holder's inequality with respect to the measure fi^idr; d£, drj). 
By virtue of (|2.22p . the linear growth of a and (jl.8p we obtain 



SUP T li3 ,2,2,2 < Ch°P, 
(t,i)e[0,T]xD 



(2.23) 



with a GlO 



2-/3 [ 
' 2 



Finaly, the estimates (I2TT5]) . (f2TT7|) and ([233]) imply that 



sup H-B^i, x) ||p < Ch a , 

(i,x)e[0,T]xD 



and a fortiori 



sup 

(t,x)e[0,T]xD 



D^u(t,x),B h Jt,x) 



< Ch c 



(2.24) 



(2.25) 



with a g]0, 71 A 72 A — ^ A ^r[. This finishes the analysis of the || • || p 
contribution to the left-hand side of (|2.1ip . 

We next consider the L p (f2, ?^r)-norm of D(D.^u(t,x),B h ^t,x)) . 

As in the previous step, we shall replace B h ^(t, x) by B h ^(t, x) with arbitrary 
t G [0,T], x G M 3 . By virtue of (HTTPl) and (12351) it suffices to study the 
L P (Q, 7-^)-norm of DB^(t,x;t,x) for i = 1,2. We start with the analysis 
of B$(t,x;t,x). 

By applying the differential rules of Malliavin calculus we have 
D.*B%(t, x; i, x) = G(t - ., x - *)</(«(.*) )£>(s, y; i, x) 



+ 



G(t-s,x-y) a (u(s,y))D.*u(s,y)V(s,y;t,x) 



+a'(u(s, y))D.*V(s, y; t, x)] Af (da, dy). 
Applying Holder's inequality and using that a' is bounded, we obtain, as in 



- .,x - *),r / (ti(.*))D(s,y;t,x)| 



H7 
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<C I ds sup {E\V(s,y;t,x)\ p } . 
Jo veKP(s) 



(2.26) 



For fixed t, x we consider the Wr-valued process defined by 



K(s, y; t, x) = a" (u(s, y))D,*u(s, y)V(s, y; t, x) + a'(u(s, y))D.*V(s, y; t, x), 

(2.27) 

(s,y) € [0,T] x M 3 , for which we have 



E 



(\\K(s, y; t, x)\fn T ) < Ik" IISo V\ I *)) 2P ) * (e \\Du(s, y) ||* 

+ || ( 7 / ||S S||DP( S ,y;i,i)||^ T . (2.28) 



We can apply the L p -estimates for stochastic integrals with respect to the 
Gaussian process M of Hilbert-valued integrands (see Equation (6.8) of 
Theorem 6.1 in [16] and [12], pg. 289) yielding 



G(t — s, x — y)K(s, y; t, x)M(ds, dy) 



LP(n,n T ) 



<C ds sup {(E(V(s, y] t,x)) 2p Y + \\DV(s, y] t,m P LPi n-nr)l- 
Jo yeK°(s) V v J ' T 7 

(2.29) 

By taking t = t and x = x and considering the inequalities (|2.26|) . (|2.29|) . 
we obtain 



\DB%(t,x^ p 



<C ds sup E {D.^u{s,y),B h Jt,x)) nT 
Jo yeKD(s) 



+ E 



{D.^u(s,y),B h Jt,x)) 



2p\ 2 



+ 



D(D.^u{s,y)iB h At,x)) HT 



Lp(Q;Ht) 

(2.30) 



Then, (|2.25j) and Gronwall's lemma yield 

D(D^u(t,x),B h Jt,x] 



sup sup 



< C sup sup 



DB?(t,x) 



LP(n;H T ) 



Lp(n-,n T ) 
+ h ap ) , 



(2.31) 



with a €]0,7i A 72 A ^ A ±±^[. 
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The last step of the proof consist of checking that for an arbitrary 
bounded set DcR 3 , 



sup 

(t,i)£[0,T]xD 



DB?(t,x) 



LP(n-,H T ) 



< Ch a , 



(2.32) 



2-/3 A 1+8 [ 



with a G]0, 71 A 72 A 2 " 2 

The proof of this fact can be done following the same lines as for (|2.24p . 
We apply the results on the densities u(dr; d£,;dr)), n±(dr; d£; dr]), /x^ (dr; d£; dr]), 
defined in (|2.14p . (|2.16p . respectively, proved so far. Instead of the process 
{a(u(s, y)), (s, y) € [0, T] x M 3 } and the L p (f2)-norm, we shall deal here 
with the 7^T - valued process {D (a(u(s, y))) , (s, y) € [0, T] x M 3 } and the 
L p (f2; 7-^r) - norm. In addition to (|1.7p . we should also apply (|1.16jl and 
(|1.14[) . We leave the details to the reader. 

Together with (|2.3p and (|2.8j) this proves (|2.2p and concludes the proof 
of the first step of the proof. 



Remark 2.2 Applying first hi. 12^ and then estimates for the \\ ■ \\2, p -norm 
of the stochastic integral leads to 



Ti 



3,2 



< c 



JR 3 



a(u(s, y)) [G(t + h - s,x — y) - G(t - s.x - y)} M(ds, dy) 



2,4 3 



Thus, we loose accuracy. This may be a justification of a pretty tricky ap- 
proach in the preceding proof. 



The rest in the time expansion 

The second and last term in (12.11) to be examined is 



R(t, x, h) 



E 



g"(u(t, x, h)) (u(t + h, x) — u(t, x)Y 



We shall apply (|1.12p to the random variables £ := u(t,x,h) and Z := 
(u(t + h, x) — u(t, x)) 2 . For this, we have to make sure that the assumptions 
of Lemma 11.31 are satisfied. For Z := (u(t + h, x) — u(t, x)) 2 , and the two 
choices of £ -u(t,x) and u(t + h, x)- this has been proved in [13) . Then it 
suffices to remark that the norm || • ||% T as well as || • define convex 
functions and use the definition of u(t, x, h) to conclude. 
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Consequently, 



R{t,x,h) < C\\u(t + h,x) -u(t, x)\\l A3 . 
Owing to (|1.7p and (|1.16p we conclude that 

sup R(t,x,h) <Ch 2a . (2.33) 

(t,x)e[0,T]xD 

with q G]0,7i A 72 A A ±±^[. 

The estimates (|2.2j) and (|2.33[) show that 

sup \E [g(u{t + h,x))- g(u{t, x))] \ < Ch a , 

(t,x)e[0,T]xD 

with a E]0, 71 A72 A A Therefore the mapping i €]0, T[i— > pt <x (y) is 
Holder continuous of degree a E]0, 71 A 72 A A uniformly in y G M 3 
varying on bounded sets. 

Step 2: space increments 

Fix t s]0, T] and consider the Taylor expansion 

E [g (u(t, x)) - 9 (t*(t, x))] = £ [</ (u(t, x)) («(t, x) - u(i, x))] 



g"(u(t, x, x)) x) — u(i, x)) z 



(2.34) 



where x,x £ D and x, x) denotes a random variable lying on the segment 
determined by u(t,x) and it(i, x). 

First order term 

Our aim is to prove that 

sup \E [g (u(t,x)) (u(t,x) - u(t,x))] I < C |x — x\ a , (2.35) 
te[o,T] 

with q G]0,7i A 72 A ^ A ±±^[. 

As for the time increments, we consider Equation (jl.6p and write 

3 

E [g' (u(t, x)) x) - u(t, x))] = ^ x, x), 



i=l 



with 



5i(t, x, x) = 2? <?' (u(t, x)) ( — G(i) * w + G(t) * v ) (x) 
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-G(t) * v + G(t) * v ) (x) 



S 2 (t,x,x) = E g' (u(t,x)) 



x ds G(t — s,dz)[b(u(s,x — z)) — b(u(s,x — z))] 



S 3 (t,x,x) = E g' (u(t,x)) 



[G(t 



s, x 



z) -G(t-s,x- z)] a(u(s, y))M(ds, dz) 



Let us consider Si(t,x,x). As for the term Ti ; i(t,x,h), we first apply 
the inequality (I1.12p and notice that 



Z(t; x, x) 



dt 



G(t) * v + G(t) *vq)(x)-[ j t G{t) * v + G(t) * v ) (x) 



is deterministic. Thus, it suffices to estimate the absolute value of the ran- 
dom variable Z(t; x, x) defined before. For this, we apply Lemmas 4.2 and 
4.4 in [5] which tell us that the fractional Sobolev norm of any integration 
degree p > 2 and differential order p < 71 A 72 is bounded. Hence, since p is 
arbitrary, by the Sobolev embedding theorem we have that 

sup Si(t, x,x) < sup C\Z(t; x, x)\ < C\x — x\ p , (2.36) 
t£[0,T] te[o,T] 

with p < 71 A 72 . 

We continue the proof with the study of the term S2(t,x,x). By virtue 
of (|1.12p . it suffices to find an upper bound of 



ds / G(t — s, dz) [b(u(s, x — z)) — b(u(s, x — z))] 



2,4 2 



in terms of a power of \x — x\. 

The measure on [0, t) x M 3 defined by ds G(t — s, dz) is finite. Hence, we 
can apply Minkowski's inequality and obtain for any p£ [1, oof 



ds / G(t — s, dz) [b(u(s, x — z)) — b(u(s, x — z))] 



2,P 

< ds G(t — s, dz) \\b(u(s, x — z)) — b(u(s, x — z)) || 2 

< C\x-x\ a , 
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with a £]0,7i A 72 A A ^^[. The last inequality is obtained by using 
that b and its derivatives are Lipschitz continuous and bounded functions, 
and by applying (jl.7p and (|1.16p . 
Hence, 

sup S 2 (t,x,x) < C\x-x\ a , (2.37) 
te[o,T] 

for any a <G]0, 71 A 72 A A 

To analyze S^t, x, x) we proceed in a similar manner as for Ti^(t, x, h) 
by applying first Proposition 3.9 in [16J and then (|1.12|) . We obtain 

S 3 (t,x,x) = E g"(u{t,x)){D^u(t,x),a(u(.,*)) 
x [G(t -.,x-*)-G(t-.,x- *)] l ]0 ,t] (■)}«, 



< C 



(A,*«(t, s), <r(«(., *)) [G(t - ., 5 - *) - G(t - .,x - *)} l m (•)) 



3,43 



Notice that the last expression has a similar structure than the right- 
hand side of (USD where B\{t, x) := G(t + h - •, x - *) - G(t --,x-*) is 
replaced by G(t — .,x — *) — G(t — ., x — *). Hence, we can proceed as in the 
analysis of the time increments to see that it suffices to deduce an estimate 
for 

\\(G(t - .,x- *)a(u(.,*)),a(u(.,*)) [G(t - .,x- *) - G(t - .,x- *)])« T || 3)P , 

for any p £ [1, 00 [. 

To pursue the proof, we split the argument of the above expression into 
two terms 



S3,l(t,X,x) 



S3,2(t,X,x) 



[ dr [ d£ [ drjG(t -r,x-£) 
Jo Jr 3 Jr 3 

x a(u(r, £))/(£ ~ r])<j{u(r, rf))G{t -r,x-rj) 

f dr f d£ I drjG(t-r,x-£,) 
Jo Jr 3 Jr 3 

x a(u(r, £))/(£ - rj)a{u(r, rj))G(t -r,x-rj), 



and we apply the change of variables (£ — ► x — £, 17 — ► x — rj), (£ — > x — £,r) 
x — 77), respectively. We obtain 

S 3 ,l(t,X,x) - S3,2(t,X,x) 



'0 



dr / / G(t-r,d£)G(t-r,d7i)f(Z-r})a(u(r,x-£)) 
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x (a(u(r, x — 7/)) — a(u(r, x — 77))) 

1 dr [ [ G(t- r, d£)G(t - r, drj) [f(x - x - (£ - 77)) - /(£ - r?)] 

x a(u(r,x - £))a(u(r,x — 77)) (2.38) 

By Minkowski's inequality the || • ||fc iP -norm of the first term in the right-hand 
side of (|2.38p is bounded by 

'dr I [ G(t-r,dS)G(t-r,dr,)f(S-r,) 



x \\<r(u(r, x - £)) [<r(u(r, x-rf))- a(u{r, x - rf))] || fc)P 
<C\x-x\ a , (2.39) 

where the very last upper bound follows from (|1.7|) . (|1.16|) and (|1.14|) . 

For the second term of the right-hand side of (|2.38[) we apply Lemma 
6.1 of [5] which implies 

sup f dr f [ G(r,d£)G(r,dr])\f(x-x-(t-r]))-f(t-r})\<C\x-xf, 



tE[0,T] JO 

with a G]0, (2 — f3) Al[. From this and the properties (jl.8|) . (|1.14p . we obtain 
the upper bound C\x — x\ a . 
Hence we conclude 

sup S 3 (t, x, x) < C\x - x\ a , (2.40) 
te[o,T] 

for any a g]0, 71 A 72 A ^ A With (pOo^ - (pl40|) . we have proved 

The rest term in the space expansion 

The contribution of the second order term in (|2.34p comes from the 
estimate 



E 



g"(u(t, x, x)) (u(t, x) — u(t, x)Y 



which is a consequence of (|1 .7|) . 

Hence, we have proved that for any fixed y € R 3 the mapping 
Pt,x(y) is Holder continuous of degree a G]0, 71 A 72 A A uniformly 
inte]0,T]. 

The proof of the theorem is now complete ■ 
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